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ABSTRACT 

The low density hydrogen and heUum in the IGM probed by QSO absorption hnes is 
sensitive to the aniphtude and spectral shape of the metagalactic UV background. We 
use realistic H i and He ii Lya forest spectra, constructed from state-of-the-art hydro- 
dynamical simulations of a ACDM Universe, to confirm the reliability of using line 
profile fitting techniques to infer the ratio of the metagalactic H i and He ii ionization 
rates. We further show that the large spatial variations and the anti-correlation with 
H I absorber density observed in the ratio of the measured He ii to H i column densi- 
ties can be explained in a model where the H i ionization rate is dominated by the com- 
bined UV emission from young star forming galaxies and QSOs and the He ii ionization 
rate is dominated by emission from QSOs only. In such a model the large fluctuations 
in the column density ratio are due to the small number of QSOs expected to con- 
tribute at any given point to the He ii ionization rate. A significant contribution to UV 
emission at the He ii photoelectric edge from hot gas in galaxies and galaxy groups 
would decrease the expected fluctuations in the column density ratio. Consequently, 
this model appears difficult to reconcile with the large increase in He ii opacity fiuctua- 
tions towards higher redshift. Our results further strengthen previous suggestions that 
observed He ii Lya forest spectra at z ~ 2 — 3.5 probe the tail end of the reionization 
of He II by QSOs. 

Key words: methods: numerical - hydrodynamics - intergalactic medium - diffuse 
radiation - quasars: absorption lines. 



1 INTRODUCTION 

As primordial density fluctuations grow through gravita- 
tional Jeans instability within CDM-like models, the web- 
like distribution of dark matter strongly infl uences the clus- 
tering of baryonic matter on large scales (|Bi_et_alJ [l992 ; 



Cenll992l:lHernauist et alllQQdiMiralda-Escude et alllQQE : 
Theuns et alJll998l: IZhang et al.lll998fl . After the epoch of 
hydrogen reionization, the hydrogen gas present in this bary- 
onic component is highly ionized, leaving only a small frac- 
tion of the gas as neutral atoms in photoionization equi- 
librium within the intergalactic medium (IGM) . This dis- 
tribution of neutral hydrogen is observable as the forest of 
Lya absorption lines in the spectra of high redshift QSOs 
(see lRauchlll99d for a review). Helium-4 is the second most 
abundant nuclide in the Univ erse after hydrogen , with a 
mass fraction of Y ~ 0.24 (e.o. lWalker et alJll99ll) . Follow- 
ing the onset of He ii reionization, the structures responsible 
for the absorption seen in the H i Lya forest should also be 
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observable as discrete He ii Lya absorption lines at the rest- 
frame wavelength of 304 A (e.g. iMiralda-Escudc fc Ostrikerl 
119901: lMiralda--Escuddll993l: iMadau fc Meiksin,.1994) . 

The development of space-based UV-capable observa- 
tories in the 1990s led to the first detection of a possible 
He II Gunn- Peterson trough in the spectrum of Q0302 — 003 
by [jakobscn ct al. ( 1994) using the Faint Object Cam- 
era on the Hubble Space Telescope (HST). Higher resolu- 
tion studies of the same spectrum using the Goddard High 
Resolution Spectrograph and the Space Telescope Imag- 
ing Spectrograph (STIS) on the HST, as well as spectral 
observations of HE 2347 — 4342 using the Far Ultravio- 
let Spectroscopic Explorer (FUSE), have all demonstrated 
a strong correlation between resolved He ii absorption fea- 
tures and those seen in the H i Lya forest, confirming the- 
oretical predictions. Absorption from He ii in the IGM has 



so far been detected in six QSO spectra; PKS 1935 — 692, 


iTvtler et al.lll995l IJakobsenlll996: Anderson et al.lll998l). 


HS 1700 + 64. iDavidsen et al. 


|l996^ HE 2347 - 4342, 


jReimers et alJ Il997l: iKriss et al 


120011 ISmette et al.ll2002l: 



IShuU et alJl2004l: IZheng et alJl2004bir SDSS J2346 



0016, 
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JZhene et al. 


l2004ari. O0302 - 003 


iHoean et alJ 


19971: iHeao et alJbOOOl 



Jakobsen et al. 



1994 
20031) 



and most recently HS 1157 - 3143, jReimers et alJl2005a^ . 
The higher redshift He ii absorption measurements [z > 
2.8) exhibit a patchy distribution, with regions of high opac- 
ity interspersed by voids. This observation, coupled with the 
steady decline in the He ii opacity at lower redshift, along 



with measurements of the IGM temperature (Theu ns et alJ 
|2Qq3)i the evolution of the H i Lya forest effective opti- 



cal depth iBernardi et al 
C iv/Si IV ratio iSongail; 



2003?) and measurements of the 
199^ . has led to the suggestion 
that the final stages of the He ii reionization epoch may have 
occurred around z ~ 3. 

In tandem with the observations, detailed models of the 
IGM have been developed to examine the He ii Lya for- 
cst from a theoretical perspecti ve (e.g. IZheng & David^gnl 
1995; Bi fc Davidsen 1997 ; Cro ft et alj lT 997: Zhang et'^ 
19971: lTheuns""et al jHtj^S^ adslcv et al. 1999; Zheng et al. 
199F ; iFardal et al.lll998l: ISokasian et al.ii2002i: .Oleser et aL 



20051) . The He ii Lya forest is of particul ar interest as it is a 
sensitive probe of the low density IGM llCroft et al.lll99'^ . 
Additionally, the ratio of He ii to H i column densities, rj = 
NB_eii/Nm, provides a constraint on the spectrum of the 
metagalactic UV background (UVB) between 13.6 eV and 
54.4 eV. Measurements of r; indicate a wide spread of values, 
from r; < 1 to > 1000 , implying that the metagalactic ra- 
diation field exhibits i nhomogeneities at the scale of ~ 1 Mp c 
jReimers et al.ll2004l : IShuU et aLlliool IZheng et alJl2004bl) . 
There is also some evidence for an anti-corre lation between 
7? and the density of t he H i Lya absorbers ||R«mer^^l 



2004'; "ShuU ct al. 2004). Semi-analvtical models llHea^^^l 
200Q: .Smette et al. 2002) suggest that a soft UVB with a sig- 
nificant stellar contribution is required to reproduce the high 
opacity regions (rj > 200), while the observed opacity gaps 
{rj < 100) have been attributed to hard sources near the line- 
of-sight creating pockets of highly ionized helium. The small- 
scale variations in rj c ould also be attributed to a spread i n 
QSO spectral indices llTelfer et alJl2002l: IScott et al 1120041) . 
local density va riations jMiralda-Esc ude et al.ll2000l) . finite 
QSO lifetimes iReimers et al. 2005b) , intrinsic absorption 
within the nucle i of active galaxies ('j akobsen et al.ll2003t 
IShull et al.]l2004^ or th e filtering of QSO radiation by ra- 
diativ e transfer effects dShull et alJl2004l : iMaseUi fc Ferraral 
[2005'). However, the exact nature of these fluctuations and 
their interpretation is unclear. 

In this paper we use realistic H i and He ii Lya for- 
est spectra, constructed from state-of-the-art hydrodynam- 
ical simulations of a ACDM Universe, to test the relia- 
bility of using line proflle fitting techniques to infer the 
UVB softness parameter, defined as the ratio of the H i and 
He II metagalactic ionization rates, S — FHi/rHcii. We ob- 
tain improved estimates of the UVB softness parameter and 
its uncertainty using published estimates of the H i and 
He II Lya forest opacity. We concentrate on the redshift 
range 2 < z < 3 where the fluctuations of the mean 
He II opacity are still moderate and He ii reionization is 
probably mostly complete. Finally, we make a detailed com- 
parison between the observed fluctuations in the column 
density ratios of correspo nding H i and He ii Lya absorp- 
tion features measured bv lZheng et al.l i2004lJ) (hereinafter 
Z04b) and those expected due to fluctuations in a meta- 
galactic He II ionization fleld dominated by emission from 



Table 1. Simulations used for our resolution and box size study. 

Name Box size Particle Mass resolution 

comoving Mpc number Mq / gas particle 



15-400 


Ibh- 


1 


2 X 400^ 


6.78 X 10^ 


15-200 


15h- 


1 


2 X 200^ 


5.42 X 10® 


30-400 


30/1- 


1 


2 X 400=^ 


5.42 X 10® 


15-100 


15h- 


1 


2 X 100^ 


4.34 X lO'^ 


30-200 


30h- 


1 


2 X 200^ 


4.34 X lO'^ 


60-400 


60/1- 


1 


2 X 400^ 


4.34 X lO'^ 


30-100 


30/1- 


1 


2 X 100=* 


3.47 X 10* 



an i2004h 



QSOs and hot gas in collapsed haloes (|Miniati et 
respectively. 

The paper is structured as follows. The method we use 
for obtaining the UVB softness parameter from our simu- 
lations is described in section 2. In section 3 we present a 
simple model for UVB fluctuations. We discuss our line pro- 
file fitting method and compare the UVB softness parameter 
inferred directly from our simulations to the value obtained 
from the column density ratio rj in section 4. In section 5 
we explore the effect of UVB fluctuations on the measured 
column density ratio. Our constraints on the UVB softness 
parameter are given in section 6, along with a discussion of 
implications for the sources which contribute to the UVB. 
We summarise and conclude in section 7. 



2 SIMULATIONS OF THE H i AND 
He II FOREST 

2.1 Numerical simulations of the IGM 

We have performed seven simulations of differing box 
size and mas s resolution using the parallel Tree SPH code 
GADGET-2 dSoringel et alJl200lHSpringe^l2005^ . We shaU 
only discuss the changes made to the simulations for this 
work; unless otherwise stated all other aspects of the code 
are as described in Bolton ot al. (2005). One major differ- 
ence is the incorporation of a non-equilibrium gas chem- 
istry solver for the radiative cooling implementation. This 
solves a set of coupled time-dependent differential equa- 
tions which govern the abundances of six species (H i, H ii, 
He I, He ii. He iii, e~) and the gas temperature fe.o. |Cen| 
1992 : lAbel et alJll997t lAnninos et all Il997t iTheuns et all 
199a) . Previously radiative cooling was implemented using 
the prescription ofp^atz ot al. ( 1996), which assumes ioniza- 
tion equilibrium. This becomes a poor approximation during 
H I and He ii reionization because gas temperatures are un- 
derestimated; less atoms are ionized per unit time compared 
to the non- equilibrium solution. In addition, the radiative 
recombination rates of Abel ct al. (1993) are adopted and 
their coUisional excitation cooling term for He i is added. 
The c oUisional ionization rates are taken from lTheuns et all 
Jl998h . 

For the UVB w e have assumed a modified version of 
iMadau et alJ il999l) model, based on contributions from 
both QSOs and galaxies. A redshift mapping of the UVB 
was made so that H i reionization begins at 2: = 6 and 
He II reionization and reheating is postponed until z ~ 3.2. 
The photoionization balance is calculated assuming the gas 
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Figure 1. Example of the artificial H I (top panel) and He II (bottom panel) Lya forest spectra used for line profile fitting. The spectra 
are constructed from the z = 2.4 output of the 15 — 200 model. Each panel is constructed from four lines-of-sight. The resolution and noise 
properties resemble those of the observed spectra of Z04b. 



is optically thin, whi ch leads to an underes t imate of the 
photo -heating rates l|Ahel fc HaehnelO Il99fll : iBolton et all 
|200J). We have thus increased the He ii photo- heating rate 
to obtain ga s temperatures consistent with observations 
jR,icotti et al.i 2QQ0: Schave et al. 20QQ: McDonald c t aL 
I2OOII) . The required increase in the He 11 photo- heating rate 
w as a factor 1.8 atz > 3.2. This is somewhat smaller than 
in lBolton et alJ (|200i'J ) due to the assumed late reionization 
of He II and the effect of the non-equilibrium solver. 



parameters consis- 
J2003D. 



We adopt cosmological 
tent with those quoted by ISpergel et alj 

(!:^m,nA,^^b/l^/l,cr8,n) = (0.26,0.74,0.024,0.72,0.90,1.0), 

and a helium mass fraction of K = 0.24. The simulations 
run for our resolution study are listed in Tabled In order to 
investigate the dependence of the softness parameter on the 
assumed cosmological model, four additional simulations 
were run with flm = [0.22,0.30] and as — [0.8,1.0] using 
200'^ gas and dark matter particles within a 15h~^ comov- 
ing Mpc box. All the other cosmological and astrophysical 
parameters adopted for these additional simulations are the 
same as for the resolution study. 



2.2 Artificial spectra and the H i and He ii 
effective optical depth 

For each simulation we have constructed artificial H i and 
He II Lya forest spectra for 1024 random lines-of-sight par- 
allel to the simulation box boundaries at eleven different 
redshifts. The spectra are generated at intervals of Az = 0.1 
within the redshift range 2 < 2 < 3. Examples of the artifi- 
cial spectra are shown in figure 1. The pixel size and signal- 
to-noise ratio are chosen to mimic the observed spectra used 
by Z04b in their analysis. The H i spectra are binned onto 
pixels of width 0.1 A and the signal-to-noise ratio is 40 
per pixel. The He 11 spectra have been binned onto pixels of 
width 0.025 A and the signal-to-noise ratio is 4 per pixel. 
Note that the pixel widths are identical in velocity space. 

As described in more detail in the last section, the sim- 
ulations have been run using the updated UVB model of 
iMadau et alJ lll999l) . As usual, we have rescaled the H i and 
He II optical depths in each pixel of the simulated spectra 
to match the respective observed effective optical depths, 
Teff = — ln(F), where {F) is the mean flux of the 1024 lines- 
of-sight. The H I and He 11 optical depths scale inversely 
with the H i and He 11 ionization rates; the metagalactic 
H I and H e 11 ionization rates can thus be determined in this 
way {e.g. iRauch et allliggTl: iMcDonald fc Miralda-Escudj 
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Figure 2. Compilation of He II Lyo forest opacity measurements 
taken from the literature. The solid curve is a fit to the data. The 
estimated uncertainty in this fit is shown by the filled circles with 
error bars. 



l200lfl . We will use this later to constrain the softness pa- 
rameter of the UVB. 

For the H i effective optical depth the central values 
from the fitting formula of Schavo ct al. ( 2003) have been 
used. The corresponding uncertainties have b een estimated 
by bin ning the Icr errors given in table 5 of ISchave et alJ 
J2003I) into redshift bins of width Az = 0.1. 

Measurements of the Ife ii effective optical depth, thcII, 
are more difficult to obtain; finding QSOs which have rel- 
atively clear lines-of-sight and which are bright enough in 
the far-UV to enable the detection of He ii absorption spec- 
tra is challenging. Only six QSO spectra are known to show 
intergalactic He ii absorption, and only four have yielded 
He II opacity measurements. It is therefore not clear to what 
extent the available data provides a statistically representa- 
tive measure of the He ii opacity. 

Nevertheless, observational studies of the He ii Lya for- 
est have met with impressive success in recent years, in par- 
ticular with high resolution studies made using STIS and 
the Far- Ultraviolet Spectroscopic Explorer (FUSE). Figure|5| 
shows a compilation of He ii optical depth measurements 
taken from the literature. The data are from Q0302 — 00 3 
Jjakobsen et al. 1994: Hoe an et alJ ll997: Heap ct al. 2000), 
PKS 1935 - 6 92 (.JakobsenI 11999: [Anderson et al.. .199al . 
HS 1700 + 64 jPavidsen et al.lll99d) and HE 2347 - 4342 
jReimers et al.lll997l: iKriss et alJl200ll: IZheng et alJl2004bl) 
with uncertainties shown where applicable. The solid curve 
in figure |5| shows the fit we use as the mean He ii effective 
optical depth for our artificial spectra, given by: 

logTHcii = 0.345 + 5.324 log ■ (1) 

Note that this fit only provides a general description for 



the global evolution of the He ii effective optical depth. The 
He II Lya forest opacity exhibits strong fluctuations which 
increase rapidly with increasing redshift; the metagalac- 
tic He II ionization rate is ob viously not s patially uniform 
jReimers et alJl2004l: IshuU et al. 2004; Zhe ng et alj|2004bll . 
The error bars attached to the solid circles in figure |21 pro- 
vide an estimate of the variation in the He ii opacity using 
a simple model for a fluctuating UVB due to emission from 
QSOs. The error bars show the 25"^ and 75**^ percentiles of 
thcII for each set of 1024 spectra at the relevant redshift. We 
will turn to the details of modelling these UVB fluctuations 
in section 13 



2.3 Numerical convergence 

As discussed by iTheuns et alj 



Jl998 h and ' Bolton et all 

12005), inferring the metagalactic ionization rate by repro- 
ducing the observed effective optical depth with artificial 
Lya forest spectra stretches present-day numerical capabil- 
ities due to the wide range of physical scales involved. It is 
thus important to perform numerical convergence tests. We 
have used the simulations listed in Table Q for a mass res- 
olution and box size study. The 400^ simulations are only 
run to z — 3; they were too computationally demanding to 
enable a practical study below this redshift. However, the 
z = 3 data should provide a good indication of how close 
the simulations are to convergence. 

As box size is increased, larger voids can be accommo- 
dated within the simulation, reducing the mean Lya ab- 
sorption and altering the H i and He ii gas distribution in 
a similar manner. There is an 8 per cent reduction in both 
the H I and He ii ionization rates inferred from the 30 — 200 
model compared to the 15 — 100 run with the same mass res- 
olution, with a further 3 per cent reduction for the 60 — 400 
data. Note that the ratio of the inferred ionization rates for 
the two species d epends little on box s ize. 

As noted by ITheuns et al.1 (Il998li . convergence of the 
He II opacity requires higher numerical resolution than the 
H I opacity. Increasing the mass resolution of a simulation 
will resolve smaller haloes, transferring more optically thin 
gas from the low density IGM into optically thick, high 
density regions, decreasing the mean Lya absorption. The 
reduction of the H i effective optical depth resulting from 
this re-distribution of the low density gas is partially offset 
by the increased contribution to the opacity from the high 
density regions. For He ii, which is optically thick down to 
much lower gas densities, this offset is less pronounced, pro- 
ducing a greater overall change in the He ii effective opti- 
cal depth for increased mass resolution. The H i ionization 
rate inferred from the 15 — 200 model is 8 per cent lower 
compared to the value of the 15 — 100 run, with a further 
drop of 8 per cent for the 15 — 400 model. The HI opac- 
ity has only marginally converged at z = 3. As expected the 
He II ionization rate shows a greater reduction in magnitude 
with increasing mass resolution, exhibiting a 20 per cent re- 
duction between the 15 — 200 and 15 -- 400 models and a 33 
per cent drop for the 30 — 200 and 30 — 400 simulations. Note 
that the stricter resolution requirement for thoII compared 
to thi has implications for the softness parameter we will 
infer later. For lower mass resolution the inferred softness 
parameter will be smaller compared to that obtained from 
identical higher resolution simulations. 
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We estimate that the ratio of the He ii and 
H I ionization rates inferred from the 15 — 200 simulation 
will be systematically low by at least 12 per cent in the red- 
shift range 2 < z < 3. We shall use this to correct our final 
estimation of S in section |S| 



3 A SIMPLE MODEL FOR UVB 
FLUCTUATIONS 

3.1 Fluctuations in the He ii opacity due to QSOs 

At 2 < 2 < 3 fluctuations in the H i metagalactic i oniza- 
tion rate are expected to be smal l {e.g. iCroft et al]|l999t 
iMeiksin fc Whitell2004l: ICroftil2004ri and the observed opac- 
ity variation is generally believed to be due to density fluc- 
tuations in the IGM. However, the spatial fluctuations in 
the He ii opacity are significantly larger than those of the 
H I opacity at the same redshift. Furthermore, the ratio 
of the He ii to H i opacity also shows large spatial fiuctu- 
ations llReimers et al.ll2004l : IShuU et alJ 120041: IZheng et all 
liooiS). The only plausible explanation for these observa- 
tions are substantial fiuctuations in the metagalactic radi- 
ation field at the He ii photoelectric edge. If, as is gener- 
ally assumed, the He ii ionizing fiux is produced by QSOs 
this conclusion is not too surprising. The comoving attenu- 
ation length for He ii ionizing photons is abo ut a factor of 
ten smaller than that for H i ioniz ing photons iFardal et alJ 
Il998l: iMiralda-Escude et alj l2000l) and thus is similar or 
shorter tha n the mean d i stance between bright QSOs. As 
discussed bv Fard al et alj lll998f) this is expected to result in 
substantial fiuctuations in the metagalactic He ii ionization 
rate. The details of these fiuctuations will depend on the 
luminosity function, life time, spectral energy and angular 
distribution of the emission of QSOs as well on details of the 
spatial distribution of He ii. A full numerical model of this 
problem including self-consistent radiative transfer around 
point sources is probably beyond p r esent-day numerical ca- 
pabilities (see ISokasian et al J 12003 : iMaselh fc Ferrarall2005l 
for attempts to implement some of the radiative transfer ef- 
fects) and certainly beyond the scope of this paper. We will 
in stead use a simple "toy model" similar to that employed 
bv lFardal et al.l 1^9^ which should incorporate many if not 
most of the relevant aspects of a fluctuating He ii ionization 
rate. 

We first assume that at any given point all QSOs within 
a sphere of radius equal to the typical attenuation length, 
-Ro, contribute to the He ii ionization rate. This assumption 
will approximate the effect of an optically thick medium on 
the propagation of ionizing radiation. Note however that a 
fully self-consistent treatment of radiative transfer will be re- 
quired to accurately model this process. The He ii ionization 
rate at redshift z can then be written as. 



rHcIl(2:) = 



Nqso 

E 



Lj{ly) (7HcIl(!^) 



hv 



dv 



(2) 



for _Rj < Ro/ (1 + z), where Rj is the proper distance of QSO j 
with luminosity Lj{v), aHeii(j') is the He ii photoionization 
cross-section and all other symbols have their usual mean- 
ing. We then populate collapsed haloes in our simulation 
with QSOs by Monte-Carlo sampling the observed QSO 



luminosity function. The comoving attenuation length of 
He II ionizing photons we use is, 

Ro = 30 {^-^) ' Mpc, (3) 

assuming the number of Lyman lim it systems per unit red- 
shift is proportional to (1 -I- zY'^ dStorrie-Lombardi et alJ 
L994j and using a normalisation based on the model of 
|Miralda-Escudc et al. (2000). Note that this length is some- 
what larger than our simulation box. We have thus tiled 
identical copies of our periodic 15-200 simulation box in all 
three spatial directions around a central "master box" to 
create a volume large enough to fully contain the "atten- 
uation sphere" around each point in the master box from 
which the artificial spectra are produced. 

We have used the fit to t he optical lu minosity function 
(OLF) of QSOs obtained by iBovle et al.) fcoOQ) in the bj 
band (effective wavelength 4580 A), 

4it/Lt,{z) 



(p{L,z) = 



[L/L.(2)]-/3i + [L/U{z 



(4) 



For the evolution of the break luminosity, L «(z) we have 
taken the pure luminosity evolution model of iMadau et alJ 
ftl99aV 



U{z) ^ L,{Q){1 + zY 



1 6*^^(1 + 6^^-) 
e?z -I- e€^* 



(5) 



The value of the model/fit parameters are C, = 2.58, ^ = 
3.16, = 1.9, /3i = -3.41, 02 = -1.58, log(V./Gpc-^) = 
2.99 and log(L. (O)/L0) = 10.67 for S^^ = 0.3, = 0.7 
and h = 0.7. A power law spectral distribution for the QSO 
luminosity, L{v) oc , is assumed. Note that more re- 
cent estimates of th e OLF exist, {e.g. ICroom et al.ll20o3 : 
iRichards et al J 120051 ), but we adopt the a bove as it is the 
basis of the updated iMadau et alJ (119991) UVB model to 
which we compare our data. 

Collapsed haloes within the simulation volume have 
been identified using a friends-of-friends halo finding al- 
gorithm with a linking length of 0.2. We have based the 
assignment of a particular QSO luminosity to a collapsed 
halo on the empirically determined relation between black 
hole mass and the veloci ty dispersion of its host bulge 
l)Ferrarese fc Merritil l200ff) . For a QSO luminosity Lqso 
randomly drawn from the OLF, the velocity dispersion of 
the host bulge is approximately: 



200 kms 



^QSO 



4.61 X 1012 e L 



(6) 



where e is the radiative efficiency of the black hole as a frac- 
tion of the Eddington limit. We adopt e = 0.1. The halo 
with the velocity dispersion closest to this value is then cho- 
sen and its position is determined by randomly selecting the 
identified halo within one of the boxes tiled to create the to- 
tal simulation volume. The total number of QSOs brighter 
than absolute magnitude Mum in our comoving simulation 
volume, V, at redshift z is then required to satisfy: 



NQSoiz,Mb, < Mlin,) = V 



0(Mb 



(7) 



where we adop t = — 22 . 

Following iMadau et alJ 1^9^, we adopt a broken 
power law for the QSO spectral energy distribution: 
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L{u) (X < u 



(2500 < A < 4600 A), 
(1050 < A < 2500 A), 
(A < 1050 A). 



(8) 



The extreme-ultraviolet (EUV) spectral index, a, is a vari- 
able parameter in our model. Each QSO placed in the simu- 
lation volume has a value for a assigned to it by Monte Carlo 
sampling the distribution of QSO spectral indic es in the 
EUV. The distribution is based on the data of Telfer et alJ 
i2002h for radio quiet QSOs, taken to be a Gaussian with 
mean 1.61 and standard deviati on 0.86. Note, h owever, 
lower redshift FUSE observations llScott et alJl20o4) favour 
a harder EUV spectral index. 

When constructing artificial spectra using our fluctuat- 
ing UVB model, the opacity in each pixel in the artificial 
He II spectra is rescaled linearly with the new fluctuating 
ionization rate computed using equation |21 The new set of 
spectra are then rescaled once more to match the observed 
He II effective optical depth. We do not alter the opacity of 
the corresponding H i spectra; a uniform UVB at the hy- 
drogen photoelectric edge should be a reasonable approxi- 
mation at the redshifts we conside r {e.g Croft ct al. 1999t 
iMeiksin fc Whitell2004l: ICrofdl200 . Note that we have also 
used this fluctuating UVB model to estimate the uncertainty 
in the He ii effective optical depth in figure |5| 



3.2 Fluctuations in the He ii opacity due to 
thermal emission from hot gas 

It has recently been suggested bv lMiniati et alJ ||2004) (M04 
hereinafter) that thermal radiation from shocked gas in col- 
lapsing haloes may provide a substantial contribution to 
the UVB at the He ii photoelectric edge. M04 show that, 
for plausible assumptions for the distribution of hot gas, 
bremsstrahlung and line emission could provide a similar 
number of He ii ionizing photons compared to QSOs at 
z — 3, and may even dominate the total photon bud- 
get above the He ii photoelectric edge for z > 4. However, 
within the average attenuation length for an He ii ionizing 
photon, there will be many more haloes contributing to the 
UVB via thermal emission compared to QSOs, which have 
a rather small duty cycle. It is therefore unclear whether 
an He ii metagalactic ionization rate dominated by emission 
from hot gas will reproduce the observed fluctuations in the 
He II opacity, or whether the more uniform distribution of 
sources expected for thermal emission will damp fluctuations 
in the metagalactic radiation field at the relevant physical 
scale. 

We modify our existing model for a fluctuating UVB 
due to QSOs to explore this. Following M04, a halo of mass 
M will have virial temperature. 



T(M) ~ 2.3 X 10^ 



M 



2/3 



(^) 



(9) 



The number of haloes emitting ionizing radiation at a given 
redshift will then be constrained by the radiative lifetime 
of the gas in each halo. Using the prescription of M04, the 
ratio of the gas cooling time to the Hubble time for a halo 
of mass Af is given by, 

^/^'1 + .X-V2 



^cool 
^Hubble 



0.21 



M 
1012 M, 



(10) 



where the last term takes into account the expected exten- 
sion to the gas cooling time due to energy injection by su- 
pernovae. We use equation 1101 to determine the probability 
that each halo within the simulation volume is radiating, 
thus providing a lower limit to the number sources within 
an attenuation sphere at any given redshift. We assume the 
monochromatic luminosity of the hot gas in the halo scales 
as MT{M)~^^'^ , and the spectrum of the bremsstrahlung 
emission is flat (/^ = constant) for the purposes of comput- 
ing FHeii- We then adjust the flux of the radiating haloes 
with a global factor such that the integrated flux repro- 
duces the observed He ii optical depth. This assumption is 
extreme in the sense that it leaves no room for the contri- 
bution by QSOs. Note, however, that our assumption for 
the mean free path of He ii ionizing photons is probably an 
underestimate for the rather hard bremsstrahlung spectrum. 



4 MEASURING THE SOFTNESS 
PARAMETER FROM H i AND 
He II Lya FOREST ABSORPTION LINES 

4.1 The softness parameter 

A combined analysis of the H i and He ii forest makes it 
possible to constrain the spectral shape of the UVB. This 
is normally done by defining the ratio of the metagalactic 
H I and He ii ionization rates as a softness parameter, 

Thi 



S = 



(11) 



This quantity provides a direct constraint on the spectral 
shape of the UVB between the H i and He ii photoelectric 
edges at 912 A and 228 A. Assuming the hydrogen and he- 
lium gas in the IGM is in ionization equilibrium, the ratio 
of observed column densities Nneii/Nm measured from ab- 
sorption lines in the H i and He ii Lya forest ca n be related 
to the softness p arameter by {e.g. iMiralda-Escud j Il993l : 
iFardal et ai]ll99d) 

A^HcII „WHcIII aHelll /.^.n 

V=^j — = S , (12) 

JVHi nHii ami 

where ni and ai are the number density and radiative re- 
combination coefficient for species i. If the IGM is highly 
ionized, 

Y 



J^HBIII 

nmi 



•n-Hc 



4(1 - Y) 



(13) 



Adopting a helium mass fraction of y = .24 and evaluating 
the radiative recombination coefficients of lAbel et alJ (|l997^ 
at T = 15000 K: 

S ~ 2.47;. (14) 

It is important to note that equation 1141 is an approxima- 
tion and as such will not provide an exact conversion from 
the observable rj to S. In particular, if He 11 is still under- 
going reionization equation 1141 will underestimate the true 
softness parameter; not all the helium in the IGM will be 
doubly ionized and equation 1131 will no longer hold. Simi- 
larly, if He II is still being reionized, the assumption of ion- 
ization equilibrium will no longer hold and equation 1141 will 
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overestimate the true softness parameter. There is also a 
weak temperature dependence in this relation depending on 
the assumed values for the recombination coefficients {e.g. 
fShuU et al.,.200 4l. However, we expect this approximation is 
accurate to within 5 per cent for the redshift range we con- 
sider. Equation 1141 will be adopted throughout this paper to 
convert from rj to S. 

4.2 Fitting line profiles to the H i and 
He II LyQ forest 

In their recent study of the helium reionization history, 
Z04b present H i and He ii Lya forest spectra of the quasar 
HE 2347 — 4342 obtained using the Very Large Telescope 
(VLT) and FUSE. They analyse the spectra by fitting Gaus- 
sian profiles to the H i lines in the VLT spectrum. Once 
the H I components are identified, counterpart He ii lines 
are defined in the FUSE spectrum, constrained to have 
the same redshift and line width as the H i lines; only 
the column density is allowed to vary. In the event that 
He II optical depth is too high to produce reasonable val- 
ues of ri, additional He ii sub-components are added with a 
line width 6hoII = 27 kms^^. The assumption of equal line 
widths is based on the analysis of a sub-sample of unblended 
H I and He ii lines by Z04b; they find a line width ratio of 
C ~ fecii/&Hi = 0.95±0.12. In this way Z04b have estimated 
the column density ratio in the redshift range 2.0 < z < 2.9. 
Z04b detect large variations in the column density ratio on 
small scales {Az ~ 0.001) and have interpreted these as spa- 
tial fiuctuations in the shape of the UVB. ShuU ot al. (2004) 
who used an independent pixel-by-pixel optical depth tech- 
nique to analyse the same observational data come to similar 
conclu sion s. However, as poin ted out bv lFechner fc ReimersI 
l|2004^ and lShull et"al] (|20o3), values of < 10 and rj > 400 
should be considered with some caution and are probably 
not real. They are likely to be caused by uncertainties intro- 
duced by background subtraction, the low signal-to-noise of 
the spectra, line saturation and blending with higher order 
Lyman series lines and metals. 

We have performed the same analysis as Z04b on our 
artificial Lya spectra. This enables us to assess if there 
is a substantial contribution from errors in the line fitting 
method to the observed variations in rj and also to determine 
how well the softness parameter can be recovered. We have 
analysed our artificial spectra using an automated version 
of the Voigt profile fitting package VPFIT (Carswell et al., 
!http://www. ast.cam.ac.uk/~rfc/vpfit. html I, modified to si- 
multaneously fit H I and He ii absorption lines. The auto- 
mated Voigt profile fitting procedure is run on 50 random 
lines of sight taken from outputs produced by the 15 — 200 
simulation. For each H i line which is identified, a corre- 
sponding He II line is fitted at the same redshift. Additional 
He II lines are added where no corresponding H i can be 
identified. Turbulent line broadening is assumed, such that 
the ratio of line widths is fixed to be ^ = &hoIi/6hi ~ 1.0, 
consistent with the Z04b result. This should be a reason- 
able assumption; Hubble broadening is expected to be the 
dominant contr ibution to the width of absorption lines in 
the Lya forest dWeinberg et al.lll998fl . Note, however, that 
additional motions in the IGM velocity field from conver- 
gent flows and shocks may produce a departure from pure 
turbulent broadening for some of the lines. The maximum 



column density allowed in the fitting procedure is fixed at 
10^^ cm~^. Similarly, the upper limit to the width of the 
fitted line profiles is set at 80 kms~^. 

4.3 Measuring the softness parameter using line 
profile fitting 

The left panel of figure |H] shows a scatter plot of the col- 
umn density ratio rj we have obtained from our artificial 
spectra at z — 2.4. A uniform background was assumed. 
When determining rj we only use line pairs which have well 
defined profile fits; any line with an error in the column 
density greater than a factor of two is rejected. The H i and 
He II column densities are tightly correlated with a moder- 
ate scatter. 

In the right panel of figure |3 we compare the softness 
parameter calculated from the column density ratio using 
eauation ll4l to the actual value used to produce the artificial 
spectra at 2 < z < 3. Again we assume a spatially uniform 
UVB. For the range of values under consideration (100 < 
S < 350), line fitting recovers the actual value of the softness 
parameter well. If we take the median column density ratio 
(open squares) the recovered value is on average 7 per cent 
less than the actual value used to produce the artificial spec- 
tra. Following Z04b, we also plot the logarithmic mean of 
the softness parameter log(S') = log(6'i)]/n (filled circles 
with error bars). The uncertainty on the logarithmic mean, 
log(A5') = [5^[log(AS'i)]^]^''^/n, is computed by propagat- 
ing the uncertainties on the H i and He ii column densities. 
The median appears to provide the better measure of the 
softness parameter; it is not weighted so heavily in favour 
of the lower column density lines which are more commonly 
identified by the fitting package. As might be expected, we 
find fitting Voigt profiles to only H i and He ii Lya absorp- 
tion features will not recover the softness parameter accu- 
rately once the He ii lines are predominantly saturated. In 
this case, only higher order He ii absorption lines or an accu- 
rate determination of the He ii effective optical depth from 
a large statistical sample of He ii absorption spectra would 
enable the recovery of S. 

Overall, we conclude that the line fitting method of 
Z04b will recover the spectral shape of the UVB along the 
line-of-sight they consider well, although one should keep 
in mind that the number of observed spectra is very small, 
and systematic uncertainties in the line fitting procedure 
means one should view the extreme values of rj w i th a certain 
amount of scepticism jFechner fc Reimer dliooi IShuU et alJ 
■2004') . It is also clear that the magnitude of the observed 
scatter in the column density ratio is much larger than our 
results for a spatially homogeneous UV background. We will 
discuss this in more detail in the next section. 



5 MODELLING THE OBSERVED 

FLUCTUATIONS IN THE COLUMN 
DENSITY RATIO 

5.1 The magnitude and physical scale of the 
fluctuations 

The column density ratio, rj = A''hcIi/A'^hi, measured by 
Z04b exhibits large variations (1 < rj < 1000) and 
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Figure 3. Left: Scatter plot of the He II against H I column density obtained by fitting Voigt profiles to artificial Lya forest spectra at 
z = 2.4. A uniform UVB has been assumed. The dashed line indicates the median column density ratio of the sample. Right: Comparison 
of the softness parameter, S, obtained by rescaling the artificial Lya forest spectra to match the observed H I and He II effective optical 
depths (solid line) with the value determined by fitting Voigt profiles to the absorption features. The filled circles with error bars show the 
logarithmic average of S and the associated uncertainty, as obtained from the column density ratio r) using eauation ll4l The open squares 
indicate the median value of S. 



the evolution of the median softness parameter is consis- 
tent with a UVB with a relatively hard spectrum, 77 < 
200. It has been suggested that spatial fluctuations in 
the UVB at the He 11 photoelectric edge could produce 
the observed He 11 opacity variations {e. g. iReimers et alJ 
I2OO4I: IShuU et aLlliooi IZheng et al.ll2004bD . So far we have 
adopted a spatially uniform UVB in the optically thin limit. 
At the H I photoelectric edge for 2 < 2 < 3 the typical 
mean free path of an ionizing photon should be much longer 
than the distance between ionizing sources; a large number 
of sources are expected to contribute to the H i ionization 
rate at any given point. Hence, as discussed previously, the 
observed variations in 77 must be caused by fluctuations in 
the UVB at the He 11 photoelectric edge. We will now use 
the model described in section|H]to investigate whether such 
fluctuations can explain the distribution in 77 observed by 
Z04b. 

The left column of figure |1| shows the column density 
ratio -q at z = [2.1, 2.4, 2.8] as measured by Z04b. The second 
column shows the same measured from our artificial spectra 
assuming a spatially uniform background. The scatter in 77 
for our model spectra is clearly significantly smaller than 
observed. The observed scatter is thus unlikely to be ex- 
plained by uncertainties in identifying the correspondence 
between H i and He 11 absorption lines or measuring their 
column densities. The third column of the plot shows 77 ob- 
tained from spectra constructed using our model for spatial 
fluctuations in the He 11 ionization rate due to QSOs. The 
scatter is significantly larger compared to the uniform UVB 
case, suggesting that UVB fluctuations do indeed play an 



important role in producing the observed variation in the 
He II opacity and softness parameter. The total scatter is 
around 2.0 dex, about 0.5 dex smaller than that of the Z04b 
data. It is not obvious how significant this rather small dif- 
ference is; Voigt profile fitting is a somewhat subjective pro- 
cedure and we have been rather strict in omitting bad fits. 
Nevertheless, it may indicate that we have over-estimated 
the attenuation length, under-estimated the scatter of spec- 
tral indices or that there is an additional contribution by 
sources which are harder and rarer than the QSOs we have 
modelled here. Additional physical effects such as radiative 
transfer, the lifetime of the sources and the opening angle for 
the ionizing emission may also play a role in increasing the 
scatter. Finally, systematic biases in the line fitting proce- 
dure of Z04b ma y artificially increase the o bserved scatter by 
around 0.5 dex iFechner fc Reimersll20oJ) . As already men- 
tioned the observed values of 77 < 10 and 77 > 400 should 
be considered wit h some caution and are probably not real 
JShuUeEaOEooi). The right column shows the result for 
our model of He 11 ionizing photon emission by hot gas. In 
this case the scatter is as small as for the uniform back- 
ground model. Such a small scatter for a UVB produced 
by hot gas is not surprising; the number of sources con- 
tributing to the He 11 ionizing UVB is more than an order 
of magnitude larger compared to the QSO emission model. 
As found for the spatially uniform UVB, the small scatter 
in the column density ratio is not consistent with the ob- 
served scatter. M04 claim that the contribution of hot gas 
to the He 11 ionizing background should increase strongly 
with increasing redshift. This result appears to be at odds 
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Figure 4. Scatter plots of the column density ratio r\ measured by Z04b (left column), from artificial spectra constructed with a uniform 
UVB at the He II photoelectric edge (second column), a fluctuating UVB due to QSOs (third column) and a UVB dominated by emission 
from hot gas in collapsing haloes (right column). The redshift ranges are indicated on each panel. The dashed line shows the median column 
density ratio of the artificial spectra produced with a uniform UVB, shown in the second column. To facilitate comparison, the same line 
is repeated in the other panels. 



with the increasing fluctuations in the He ii opacity towards 
higher redshift, suggesting that emission by hot gas actually 
contributes little to the He ii ionization rate and that the 
estimates of the bremsstrahlung volume emissivity by M04 
may have been too large. This is perhaps also not too sur- 
prising as the assumed density profile of the emitting gas 
in dark matter haloes is very uncertain. More quantitative 



statements will require detailed modelling of the UVB using 
radiative transfer which is beyond the scope of this paper. 

As an additional check on our assumptions for the 
H I and He ii effective optical depth, we plot as a dashed 
Une in figure |1] the median r\ obtained from the artificial 
spectra with a uniform UVB (second panel). The same line 
is repeated in the panels for the artificial spectra constructed 
using a fiuctuating UVB and the observations. Although the 
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Figure 5. Scatter plots of the column density ratio against rcdshift for 2.3 < z < 2.5. The observational data of Z04b are represented 
by open diamonds in both panels. In the left panel, the column density ratios determined from our artificial spectra, generated using a 
uniform UVB at 2 = 2.4, are indicated by stars. The redshifts for the simulated spectra have been shifted such that they cover the same 
redshift range as the observational data while maintaining the physical separation between data points. The actual simulation output is at 
at z = 2.4. The stars in the right panel show the column density ratios obtained from artificial spectra generated with a fluctuating UVB 
due to QSOs. 



scatter in the column density ratio is much smaller for the 
artificial spectra with the uniform UVB, the median value 
for Tj is consistent with that observed. This is encouraging, as 
it suggests that our fit for the He ii effective optical depth at 
2 < z < 3 provides a reasonable representation of the true 
He II Lya forest opacity. Note the increase in the median 
column density ratio with redshift; a softer UVB is required 
to reproduce the observed H i and He ii Lya forest opacity 
at higher redshift. 

Figure 1^ shows how rj varies along the line-of-sight and 
gives an impression of the physical scales over which the 
fluctuations are correlated. The Z04b data in the range 
2.3 < z < 2.5 are plotted with open diamonds in both pan- 
els. For comparison, results for our simulated spectra are 
shown as stars. The left panel shows the simulated data as- 
suming a spatially uniform background while the right panel 
is for our model of the UVB fluctuations due to QSOs. The 
redshifts of the simulated spectra have been shifted to cover 
the same redshift range as the observational data. How- 
ever, the actual physical distance between each data point 
for the simulated spectra has been preserved to enable a 
fair comparison between the scales over which varies. For 
the simulated spectra with a uniform UVB the data points 
strongly cluster around rj ~ 70, in strong disagreement with 
the distribution for the observational data. For the simu- 
lated spectra created with the fluctuating UVB due to QSOs 
the spread in r] reproduces the observational result remark- 
ably well, although a full treatment of radiative transfer ef- 
fects may be required t o produce UVB fluctuations on the 
smallest physical scales dShull et all20oJ : lMaselli fc Ferraral 



I2OO5I) . Note again that it is not clear whether the underes- 
timate of the number of points which see a very hard spec- 
trum, < 10 in the simulations is real. 



5.2 The anti-correlation of the column density 
ratio with H i density 

The studies of UVB flu ctuations by iReimers et alJ ||2004) 
and IShuU et al.l l)200'j) have also uncovered an anti- 
correlation between the column density ratio rj and the den- 
sity of the H I Lya absorbers. Regions of high opacity within 
the H I Lya forest, corresponding to the overdense knots, 
filaments and sheets which form as a natural consequence 
of structure formation, tend on average to see a somewhat 
harder radiation field (7; < 100). In contrast, the voids in 
the H I forest account for the majority of the larger values 
of rj, corresponding to a softer metagalact ic radiation field 
{rj > 100). In particular. IshuU et all j2004^ find up to 80% 
of the void (thi < 0.05) path length has rj > 100, as opposed 
to around 40% for higher density regions (thi > 0.05). 

Figure |S| shows a plot of the cumulative probability 
distributio n function (CPDF) for rj at z = 2.4, similar to 
figure 5 in IShuU etabl ^2004 ). A fluctuating UVB at the 
He II photoelectric edge due to QSOs was assumed. The 
solid curve is the CPDF for all identified line pairs where 
Thi > 0.05. The dotted curve shows the CPDF for under- 
dense regions, rm < 0.05. There is very good agreement 
with figure 5 in IShuU et alJ (|2004) : the underdense regions 
clearly see a softer radiation field on average. In our model 
the origin of this anti-correlation can be easily understood 
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Figure 6. The cumulative probability distribution func- 
tion (CPDF) for r; at 2 = 2.4; a fluctuating UVB at the 
He II photoelectric edge due to QSOs has been assumed. The solid 
curve shows the CPDF for H I /He II line pairs with rjji > 0.05. 
The dotted curve shows the CPDF for underdense regions (rjji < 
0.05). Low density regions in the IGM appear to see a softer radi- 
ation fifildjConsistent^ observational data (e.o. .Rcimcrs ct alj 
l2004HShull et al.ll2004) . 



from the spatial distribution of the QSOs which are respon- 
sible for the radiation at 4 Ryd. The underdense regions 
are typically further away from the massive haloes which 
host these QSOs. The flux at 4 Ryd is therefore on aver- 
age smaller in underdense regions and thus the value of r/ 
is larger. Note, however, that systematic uncertainties from 
the line fitting method w ill also cont ribute to the the ob- 
served anti-correlation OFechner fc Reimera»20o3) . 



6 THE SPECTRAL SHAPE OF THE UV 
BACKGROUND 

6.1 The average softness parameter and its 
evolution with redshift 

We obtain an estimate of the average softness parameter 
from our simulations with a spatially uniform UVB by ad- 
justing the H I and He ii ionization rates such that the artifi- 
cial spectra reproduce the observed H i and He ii Lya forest 
opacity. The result is shown as the solid curve in figure |7| 
corrected for numerical resolution, with the systematic un- 
certainties shown by the grey region. We also show the same 
result estimated from our artificial spectra constructed as- 
suming a fluctuating UVB at the He ii photoelectric edge 
due to QSO emission, shown as the dot-dashed line in fig- 
ure |7| There is some mild evidence for an increase of the 
softness parameter with increasing redshift. We now pro- 



ceed to discuss this measurement and its uncertainties in 
more detail. 

6.2 Systematic uncertainties in the estimate of 
the average softness parameter 

The value of S we infer from our artificial spectra will 
depend on our assumptions for the various cosmologi- 
cal and astrophysical parameters which alter the opac- 
ity of the Lya forest, as well as the observational val- 
ues for the H i and He ii effective optical depths we adopt. 
However, as we are interested in the ratio of the H i to 
He II metagalactic ionization rates, uncertainties in param- 
eters on which the H i and He ii opacity have the same de- 
pendence will not contribute significantly to the overall un- 
certainty in S. In this instance, uncertainties from assump- 
tions for Qhh^ and the temperature of the IGM can be 
ignored; Qhh^ determines the total baryonic content, and 
the H I and He ii recombination coefficients have the same 
temperature dependence. Consequently, altering either of 
these parameters will have the same effect on the H i and 
He II ionization rates we infer. 

In contrast, the matter density as fraction of the criti- 
cal density, Qm, and the mass fiuctuation amplitude within 
8^~^ Mpc spheres, erg, both alter the rms fluctuation am- 
plitude at the Jean s scales, changing the gas distribution 
iBolton et al.ll200^ . As already discussed. He ii is optically 
thick down to lower densities compared to H i. For larger 
Jim or (78, the He ii opacity decreases more relative to the 
H I opacity as absorbing gas is transferred to higher den- 
sity regions. The opposite is also true for a smaller value 
of Qm or as. To estimate the uncertainty on the value of S 
in the redshift range 2 < 2 < 3 from these parameters, we 
adopt fiducial values of f2m = 0.26±0.04 and erg — 0.9±0.1, 
and run four extra simulations to determine S using param- 
eter values at the upper and lower end of the uncertainties. 

Assumptions for the slope 7 of the pow er-law effective 
equat ion of state for the low density gas jHui fc GnedinI 
Il997^ will also alter S. For a flatter effective equation of 
state, a larger fraction of the optically thick He 11 gas will 
have a higher temperature compared to the H i gas, reduc- 
ing the inferred He 11 ionization rate to a greater extent than 
the H I ionization rate, increasing the inferred softness pa- 
rameter. The converse is also true; a larger value of 7 will 
increase the He 11 opacity substantially more than that for 
H I, lowering the inferred S. The plausible range for 7 is 
1.0 < 7 < 1.6 faui fc GnediJl997l : lschave et al.ll200ol) . We 
estimate the effect of 7 on 5 by artificially rescaling the 
effective equation of state of 15 — 200 simulation data by 
pivoting the temperature-density relation around the mean 
gas density. 

Finally, the uncertainty in the adopted values for thi 
and thoII must be considered. We estimate th e uncertainty 
in Thi by binning the la measurement errors of lSchave et alJ 
(l2003h as described in section 2.3. We use the error bars 
shown in figure |5| to estimate the uncertainty in our as- 
sumption for the He 11 effective optical depth. These have 
been computed using the fiuctuating UVB model due to 
QSOs described in section|3 the lower and upper error bars 
correspond to the 25"" and 75"" percentiles of thcII taken 
from each set of 1024 spectra at the relevant redshift. An 
additional uncertainty of 10 per cent is added to take into 
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Table 2. Percentage error budget for the softness parameter 
at 2 = [2.1,2.4,2.8] from estimates of various cosmological and 
astrophysical parameters, listed approximately in order of impor- 
tance. The total error is obtained by adding the individual errors 
in quadrature. Our estimate for the softness parameter at these 
redshifts and the corresponding uncertainty is listed in the last 
row. A spatially uniform UVB has been assumed. 



Parameter 


z = 2.1 


z = 2.4 


z = 2.8 


THell 


~±^\ per cent 


per cent 


^4^* per cent 


Thi 


^18 per cent 


per cent 


_lg per cent 


7 


lot; 

per cent 


per cent 


per cent 


Numerical 


±10 per cent 


±10 per cent 


±10 per cent 




per cent 


±9 per cent 


±10 per cent 




per cent 


±6 per cent 


per cent 


Total 


ilg per cent 


^49 per cent 


-50^ per cent 


S 


139+f, 


19611^" 


o„,+576 
•^'^^-151 



account the marginal convergence of the simulations. The 
final error budget for S at 2 = [2.1, 2.4, 2.8] is summarised 
in Table |5| The dominant uncertainty is the measurement 
of the He 11 effective optical depth. 

6.3 Comparison to the UVB model of Madau et 
al. and other observational estimates 

In figure Q we compare the results of our analysis of the 
average softness parameter inferred from observed H i and 
He II effective optical depths with the UVB model of Madau 
et al. (1999) (hereinafter MHR99) and other observational 
estimates of 5*. The dotted and dashed curves show the 
soft ness parameter of t he updated UVB models of MHR99 
(see iBolton et aflbOOSl for a detailed description) for QSOs 
alone and QSOs plus young star forming galaxies (YS- 
FG s). Typically S < 100 for QSO dom inated UVB mod- 
els jHaardt fc Madaulll99d iFardal et al.llT99ai) . although if 
there is a wide range in QSO spectral in dices a QSO dom i- 
nated UVB could be as soft as S' ~ 400 jTelfer et al.ll2002l) . 
Due to their softer spectra, YSFGs only provide a signifi- 
cant boost to the UVB amplitude at the H i photoelectric 
edge, so that typically one expects S > 200 for a UVB with 
contributions from QSOs and YSFGs. 

Our results for the spatially uniform UVB model agree 
very well with the softness parameter of the updated MHR99 
model with contributions from QSOs and YSFGs. At the 
lowest redshift z ~ 2 there is some indication that the 
MHR99 model may slightly overpredict the softness param- 
eter. Note that the MHR99 UVB model also assumes a 
spatially uniform background, allowing a fair comparison. 
Our analysis thus adds to the growing evidence that the 
amplitude of the UVB at the H i photoelectric edge can- 
not be accounted for by QSO emissivities alone; a sub- 
stantial contribution from YSFGs seems to be required 





Giroux & Shulllll997l: IZhane et alJll997l: iTheuns et al. 


199S 


iBianchi et alJl200ll: iKriss et al.ll200ll: iHaehnelt et al 


"2001' 


"Smotto ot al. '2002'; 'Zhena- ct al." ■2004b'; 'ShuU ct al 


[2004 


Aguirre ct al. 2004; Bolton ct al. 2005: Kirkman ct al. 


I2OO5I 


. Note, however, that as discussed before, UV fluctua- 



tions will play a significant role in altering S on local scales, 
especially at z > 2.8. Some regions of low opacity are obvi- 
ously subject to a harder ionizing radiation field, producing 
values of 5* outside our error range. Similarly, overdense re- 
gions which are optically thick to He 11 ionizing photons can 
be subject to a softer radiation field. 

The spatial fiuctuations also affect the average 
He II ionization rate and thus the inferred average softness 
parameter. The dot-dashed curve in figure |7| shows the av- 
erage softness parameter obtained from our model of a spa- 
tially fluctuating UVB due to QSO emission. The aver- 
age softness parameter inferred for the spatially fluctuat- 
ing He II ionization rate is systematically smaller (0.15 dex) 
than that obtained for a spatially uniform UVB. This is be- 
cause in regions close to QSOs, corresponding to voids in 
the He 11 Lya forest, raising the He 11 ionization rate fur- 
ther will not decrease the effective optical depth to a great 
extent once the local optical depth is already low. The 
average He 11 ionization rate will thus be biased low if a 
uniform background is assumed. This effect is also found 
for UVB fi uctuations at the hydrogen ionization edge {e.g. 
iGnedin fc H amihon 2002; M eiksin fc White 2004). 

We also plot the data of Z04b and iHeap et al.l ll2000l) 
in figure [7| The logarithmic mean of the softness parame- 
ter calculated from the Z04b column density measurements 
is shown by the filled circles with error bars. The median 
value in each bin is indicated with an open square. We have 
only used data where both the H i and He 11 components 
have measured column densities; the analysis of our artifi- 
cial spectra using Voigt profile fitting indicates that these 
measurements should recover the UVB spectral shape well 
for the line-of-sight considered, with the median providing 
a better estimate of S. The results of Z04b agree very well 
with our analysis for a fiuctuating He 11 ionization rate. Note 
that the results of Z04b cannot be compared directly to the 
predictions of the MHR99 model which assumes a spatially 
uniform UVB. 



7 CONCLUSIONS 

We have used state-of the art hydrodynamical simulations to 
obtain improved measurements of the spectral shape of the 
ionizing UV background from the H i and He 11 Lya forest 
and to investigate the origin of the large spatial fluctuations 
observed in the He 11 to H i column density ratio. 

Using artificial absorption spectra we have shown that 
the softness parameter of the UV background can be accu- 
rately inferred from the column density ratio of He 11 and 
H I absorption lines obtained by applying standard Voigt 
profile fitting routines. Uncertainties in the identification 
of corresponding H i and He 11 absorption and errors in the 
determination of the column densities contribute little to 
the large fluctuations observed in the He 11 and H i column 
density ratio, although one should be cautious of extreme 
values of r;. These fluctuations must be due to genuine spa- 
tial variations in the spectral shape of the UVB. A model 
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Figure 7. Our determination of the average softness parameter of the UVB compared to the model of MHR99 and other observational 
estimates. The solid curve shows our results S corrected for resolution effects assuming a spatially uniform UVB. The grey shaded area 
indicates the theoretical uncertainty. The dot-dashed curves shows our result for the case of a spatially fluctuating He II ionization rate. The 
dotted and dashed curves represent the value of S predicted by the updated UVB models of MHR99, based respectively on contributions 
from QSOs alone and QSOs plus YSFGs. The filled circles with error bars show the results of Z04b calculated from the logarithmic mean 
of the measured column density ratios rj using eQuation ll4l The open squares show the median value of S in each redshift bin for the same 
data. Finally, the filled triangle represents the determination of S at ^ ~ 3 by Heap et al. (2000). 

where the He ii ionization rate fluctuates due to variation 
in the number, luminosity and spectral shape of a small 
number of QSOs reproduces the observed spatial variations 
of the He ii and H i column density ratio, as well as the 
observed anti-correlation of the column density ratio with 
the H I density. This is in good agreement with previous 
suggestions that the observed He ii Lya forest spectra at 
2 ~ 2 — 3.5 probes the tail end of the reionization of He ii by 
QSOs. The large fluctuations observed in the He ii and 
H I column density ratio argue strongly against a significant 
contribution of emission by hot gas to the He ii ionization 
rate at 2 < z < 3. 

We further obtain new constraints and error estimates 
for the mean softness parameter of the UVB and its evo- 
lution using the observed effective optical depth of the 
He II and H i Lya forest. We find that the inferred soft- 
ness parameter is consistent with a UVB with contributions 
from QSOs and YSFGs within our uncertainties. The dom- 
inant contribution to the uncertainty on S comes from the 
He II effective optical depth. Our results for the mean soft- 
ness parameter add to the growing body of evidence suggest- 
ing that a UVB dominated by emission from QSOs alone 
cannot reproduce the observed effective optical depth of the 
Lya forest in QSO spectra. 



ACKNOWLEDGEMENTS 

We thank Volker Springel for his advice and for making 
GADGET-2 available. We are also grateful to Francesco 



Haardt for making his updated UV background model avail- 
able to us. JSB thanks John Regan for assistance with run- 
ning simulations, and Natasha Maddox and Paul Hewett for 
helpful discussions. This research was conducted in cooper- 
ation with SGI/Intel utilising the Altix 3700 supercomputer 
COSMOS at the Department of Applied Mathematics and 
Theoretical Physics in Cambridge. COSMOS is a UK-CCC 
facility which is supported by HEFCE and PPARC. We 
also acknowledge support from the European Community 
Research and Training Network 'The Physics of the Inter- 
galactic Medium'. JSB, MGH, MV and RFC thank PPARC 
for financial support. We also thank the referee, Mike ShuU, 
for a very detailed and helpful report. 



REFERENCES 

Abel T., Anninos P., Zhang Y., Norman M. L., 1997, 

NewA, 2, 181 
Abel T., Haehnelt M. G., 1999, ApJ, 520, L13 
Aguirre A., Schaye J., Kim T., Theuns T., Rauch M., Sar- 
gent W. L. W., 2004, ApJ, 602, 38 
Anderson S. F., Hogan C. J., Williams B. F., Carswell 

R. F., 1999, AJ, 117, 56 
Anninos P., Zhang Y., Abel T., Norman M. L., 1997, 

NewA, 2, 209 
Bernard! M., et al. 2003, AJ, 125, 32 
Bi H., Davidsen A. F., 1997, ApJ, 479, 523 
Bi H. G., Boerner G., Chu Y., 1992, A&A, 266, 1 
Bianchi S., Cristiani S., Kim T.-S., 2001, A&A, 376, 1 



© 0000 RAS, MNRAS 000, 000-000 



14 J.S. Bolton et al. 



Bolton J., Meiksin A., White M., 2004, MNRAS, 348, L43 
Bolton J. S., Haehnelt M. G., Viel M., Springel V., 2005, 

MNRAS, 357, 1178 
Boyle B. J., Shanks T., Groom S. M., Smith R. J., Miller 

L., Loaring N., Heymans G., 2000, MNRAS, 317, 1014 
Gen R., 1992, ApJS, 78, 341 
Groft R. A. G., 2004, ApJ, 610, 642 

Groft R. A. G., Weinberg D. H., Katz N., Hernquist L., 

1997, ApJ, 488, 532 
Croft R. A. G., Weinberg D. H., Pettini M., Hernquist L., 

Katz N., 1999, ApJ, 520, 1 
Groom S. M., Smith R. J., Boyle B. J., Shanks T., Miller 

L., Outram P. J., Loaring N. S., 2004, MNRAS, 349, 1397 
Davidsen A. P., Kriss G. A., Wei Z., 1996, Nature, 380, 47 
Fardal M. A., Giroux M. L., ShuU J. M., 1998, AJ, 115, 

2206 

Fechner G., Reimers D., 2004, preprint ( astro-ph/0410622 ) 

Ferrarese L., Merritt D., 2000, ApJ, 539, L9 

Giroux M. L., ShuU J. M., 1997, AJ, 113, 1505 

Gleser L., Nusser A., Benson A. J., Ohno H., Sugiyama N., 

2005, MNRAS, 361, 1399 
Gnedin N. Y., Hamihon A. J. S., 2002, MNRAS, 334, 107 
Haardt P., Madau P., 1996, ApJ, 461, 20 
Haehnelt M. G., Madau P., Kudritzki R., Haardt F., 2001, 

ApJ, 549, L151 
Heap S. R., Williger G. M., Smette A., Hubeny I., Sahu 

M. S., Jenkins E. B., Tripp T. M., Winkler J. N., 2000, 

ApJ, 534, 69 

Hernquist L., Katz N., Weinberg D. H., Miralda-Escude J., 

1996, ApJ, 457, L51 
Hogan C. J., Anderson S. F., Rugers M. H., 1997, AJ, 113, 

1495 

Hui L., Gnedin N. Y., 1997, MNRAS, 292, 27 
Jakobsen P., 1996, in Benvenuti P., Macchetto F.D., 
Schreier E.J., eds. Science with the Hubble Space Tele- 
scope - II. Space Telescope Science Institute, Baltimore, 
p.153 

Jakobsen P., Boksenberg A., Deharveng J. M., Greenfield 
P., Jedrzejewski R., Paresce F., 1994, Nature, 370, 35 

Jakobsen P., Jansen R. A., Wagner S., Reimers D., 2003, 
A&A, 397, 891 

Katz N., Weinberg D. H., Hernquist L., 1996, ApJS, 105, 
19 

Kirkman D., et al. 2005, MNRAS, 360, 1373 

Kriss G. A., et al. 2001, Science, 293, 1112 

Madau P., Haardt F., Rees M. J., 1999, ApJ, 514, 648 

Madau P., Meiksin A., 1994, ApJ, 433, L53 

Maselh A., Ferrara A., 2005, preprint (astro-ph/0510258 I 

McDonald P., Miralda-Escude J., 2001, ApJ, 549, Lll 

McDonald P., Miralda-Escude J., Ranch M., Sargent 

W. L. W., Barlow T. A., Gen R., 2001, ApJ, 562, 52 
Meiksin A., White M., 2004, MNRAS, 350, 1107 
Miniati F., Ferrara A., White S. D. M., Bianchi S., 2004, 

MNRAS, 348, 964 (M04) 
Miralda-Escude J., 1993, MNRAS, 262, 273 
Miralda-Escude J., Gen R., Ostriker J. P., Ranch M., 1996, 

ApJ, 471, 582 

Miralda-Escude J., Haehnelt M., Rees M. J., 2000, ApJ, 
530, 1 

Miralda-Escude J., Ostriker J. P., 1990, ApJ, 350, 1 
Ranch M., 1998, ARA&A, 36, 267 
Ranch M., et al. 1997, ApJ, 489, 7 



Reimers D., Fechner G., Hagen H.-J., Jakobsen P., Tytler 

D. , Kirkman D., 2005a, A&A, 442, 63 

Reimers D., Fechner G., Kriss G., ShuU M., Baade 
R., Moos W., Songaila A., Simcoe R., 2004, preprint 
(astro-ph/0410588) 

Reimers D., Hagen H.-J., Schramm J., Kriss G. A., ShuU 
J. M., 2005b, A&A, 436, 465 

Reimers D., Kohler S., Wisotzki L., Groote D., Rodriguez- 
Pascual P., Wamsteker W., 1997, A&A, 327, 890 

Richards G. T., et al. 2005, MNRAS, 360, 839 

Ricotti M., Gnedin N. Y., ShuU J. M., 2000, ApJ, 534, 41 

Schaye J., Aguirre A., Kim T., Theuns T., Ranch M., Sar- 
gent W. L. W., 2003, ApJ, 596, 768 

Schaye J., Theuns T., Ranch M., Efstathiou G., Sargent 
W. L. W., 2000, MNRAS, 318, 817 

Scott J. E., Kriss G. A., Brotherton M., Green R. F., 
Hutchings J., ShuU J. M., Zheng W., 2004, ApJ, 615, 135 

ShuU J. M., Tumlinson J., Giroux M. L., Kriss G. A., 
Reimers D., 2004, ApJ, 600, 570 

Smette A., Heap S. R., Williger G. M., Tripp T. M., Jenkins 

E. B., Songaila A., 2002, ApJ, 564, 542 

Sokasian A., Abel T., Hernquist L., 2002, MNRAS, 332, 
601 

Songaila A., 1998, AJ, 115, 2184 
Spergel D. N., et al. 2003, ApJS, 148, 175 
Springel V., 2005, preprint (astro-ph/0505010l 
Springel V., Yoshida N., White S. D. M., 2001, New A, 6, 
79 

Storrie-Lombardi L. J., McMahon R. G., Irwin M. J., Haz- 
ard C., 1994, ApJ, 427, L13 

Teller R. G., Zheng W., Kriss G. A., Davidsen A. F., 2002, 
ApJ, 565, 773 

Theuns T., Leonard A., Efstathiou G., Pearce F. R., 

Thomas P. A., 1998, MNRAS, 301, 478 
Theuns T., Schaye J., Zaroubi S., Kim T., Tzanavaris P., 

CarsweU B., 2002, ApJ, 567, L103 
Tytler D., Fan X.-M., Buries S., GottreU L., Davis C., 

Kirkman D., Zuo L., 1995, in Meylan G., ed, Proc. of 

ESO Workshop, QSO Absorption Lines. Springer- Verlag, 

Berlin, p. 289 

Wadsley J. W., Hogan G. J., Anderson S. F., 1999, in HoU 
S., Smith E., eds, AIP Gonf. Proc. 470: After the Dark 
Ages: When Galaxies were Young (the Universe at 2 < z 
< 5). American Institute of Physics Press, p. 273 

Walker T. P., Steigman G., Kang H., Schramm D. M., Olive 
K. A., 1991, ApJ, 376, 51 

Weinberg D. H., Hernquist L., Katz N., Groft R., Miralda- 
Escude J., 1998, in Petitjean P., Gharlot S., eds, Proc. 
13th lAP CoUoq., Structure and Evolution in the Inter- 
galactic Medium. Editions Frontieres, Paris, p. 133 

Zhang Y., Anninos P., Norman M. L., Meiksin A., 1997, 
ApJ, 485, 496 

Zhang Y., Meiksin A., Anninos P., Norman M. L., 1998, 
ApJ, 495, 63 

Zheng W., Ghiu K., Anderson S. F., Schneider D. P., Hogan 
C. J., York D. G., Buries S., Brinkmann J., 2004a, A J, 
127, 656 

Zheng W., Davidsen A., 1995, ApJ, 440, L53 

Zheng W., Davidsen A. F., Kriss G. A., 1998, AJ, 115, 391 

Zheng W., et al. 2004b, ApJ, 605, 631 (Z04b) 



© 0000 RAS, MNRAS 000, 000-000 



